Gate-controlled one-dimensional channel on the topological surface 
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We investigate the formation of the one-dimensional channels on the topological surface under 
the gate electrode. The energy dispersion of these channels is almost linear in the momentum with 
the velocity sensitively depending on the strength of the gate voltage. The energy is also restricted 
to be positive or negative depending on the strength of the gate voltage. Consequently, the local 
density of states near the gated region has an asymmetric structure with respect to zero energy. 
In the presence of the electron-electron interaction, the correlation effect can be tuned by the gate 
voltage. We also suggest a tunneling experiment to verify the presence of these bound states. 



CD 



o 



> 

(N 



O 
O 



X 



PACS numbers: 73.43.Nq, 72.25.Dc, 85.75.-d 

Recent discovery of the two-dimensional (2D) quan- 
tum spin Hall systemjlj-Q, and its three-dimensional 
(3D) generalization, dubbed topological insulator, 
|l2| has established a new state of matter in the time- 
reversal symmetric systems. The topological order in 
the bulk with the gap guarantees the presence of the 
one-dimensional (ID) channels along the edge of the 2D 
sample, or the 2D metal on the surface of the 3D sample. 
These edge and surface states are protected by the time- 
reversal symmetry and the topology of the bulk gap, and 
are robust against the disorder scattering and electron- 
electron interactions. 

In topological surface state on 3D topological insula- 
tor, the electrons obey the 2D Dirac equations. This has 
been beautifully demonstrated by th e sp in- and angle- 
resolved photoemission spectroscopy [ij, [13] • To date, 
various interesting properties of the topological insula- 
tor have been predicted, particularly those relevant to 
magnetism 15 -3| and electron-electron interaction[23- 
25|. Since the surface state of topological insulator has 
a suppressed backward scattering for nonmagnetic im- 
purities, this is a promising material for designing the 
quantum curcuit. Therefore, electric control of transport 
on the surface of topological insulator is an important 
issue. In fact, the chiral edge channels are predicted to 
appear at the interface of two ferromagnets with mag- 
netization along z and —z directions, i.e., perpendicular 
to the surface. [2Q] Also, scanning tunneling microscopy 
experiments and related theories have revealed the elec- 
tronic states due to the impurity scattering or the terrace 
edge on the surface of topological insulator 27-3^. 



In this paper, we study the ID states that are formed 
on the topological surface under the gate electrode. The 
energy dispersion of these states is almost linear in the 
momentum with the velocity sensitively depending on the 
strength of the gate voltage. The energy is also restricted 
to be positive or negative depending on the strength of 
the gate voltage. Consequently, the local density of states 



near the gated region has an asymmetric structure with 
respect to zero energy. With the electron-electron in- 
teraction, the correlation effect can be tuned by the gate 
voltage. We also suggest a tunneling experiment to verify 
the presence of these bound states. 
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FIG. 1: (Color online) (Upper) schematics of the proposed 
setup. The gate electrode is attached on the surface of the 
topological insulator where ID helical mode is generated. 
A magnetic impurity placed on the gated region produces 
the backward scattering. (Lower) dispersions of surface and 
bound states. 

We consider 2D topological surface on which the ID 
gate electrode is attached as shown in Fig. [TJ The gate 
electrode is attached in the region < x < L. The inter- 
face between gated and ungated regions is parallel to the 
j/-axis. The electrons on the topological surface obey the 
2D Dirac Hamiltonian H = vpk ■ cr and under the gate. 
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FIG. 2: (Color online) The dispersion of the bound states 
as a function of ky for various U/Ef with kpL = 0.1 and 
Ef ~ VFkp- 



we add the potential U, namely H = vpk' ■ cr + C7. [33| 
Here, vp is the Fermi velocity, k{k') is a wavevector and 
(T is a vector of Pauli matrices. The dispersion relation 
is = ±VF^kl + fc2 = ±VF^k'^ + fc2 + U. Note that 
due to the translational symmetry along the y-axis, ky 
is conserved. Now, let us regard '^p' as a gap for a 
fixed ky and consider the bound states formed inside this 
gap. By setting k^ — in{n > 0) outside the gated re- 
gion to search for the bound states near the gate, with 
the above Hamiltonian, we have wave functions in each 
region represented as 



tp{x < 0) = ae 
ipiO < X < L) ^ be 

'ip{x > L) — de 
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In the above expressions, the common factor e**^"^ is 
omitted for simplicity. 

By matching the wavefunctions at the interfaces x = 
and L, we obtain 



By solving the above equation numerically, we plot the 
dispersion of the bound states as a function of ky for var- 
ious U/Ef with kpL = 0.1 and Ef = vpkp in Fig. [2] 
As seen, the bound states are generated with an almost 
linear dispersion in the momentum. The sign and the ve- 
locity of this dispersion strongly depend on the potential 
parameter U. 

In order to make the following discussion more trans- 
parent, let us take the limit [/ — > oo and L — > 
while keeping Z = UL = const. Then, we have 
vfh = —EtaiiZ(> 0), and hence E — ivFicos Z)ky and 
K = =F(sinZ)fcj,. When the sign of tanZ is fixed, then 
that of E is determined. There are two branches of the 
bound states, which constitute a ID helical mode. This 
is a physical realization of the Tomonaga model which 
is characterized by two dispersions defined for positive 
energy only (and no dispersion for negative energy). [sj] 

How does this helical mode manifest itself in observable 
quantities? In the following, we show its salient feature 
in local density of states and the current along the ID 
channel. First, we study the local density of states near 
the gated region {x < 0). The wavefunction of the bound 
states is given by 



i^six < 0) 



1 ± sin Z 



cosZ 
i{±l + smZ) 



(5) 



Note that the spins of the helical mode are opposite to 
each other. Namely, the expectation value of the Pauli 
matrices is given by (cr) || (0, ± cosZ, =p sin Z)*. 

The total local density of states can be calculated as 



p{E,x) = ^ J dk,dky\i^s\^KE^Es) 
i- j dky\i,s\^5{E-EB) 
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p^°^ +5p{E,x)+ pb{E,x) 



(6) 
(7) 



(4) with p^o) = 



and 



5p{E, x) = - -4^ de sin2 Wm ( e-^'^^ '^^'-^ J'""^ , , (8) 

' ' TT^vj, Jo \ cosOcosZ -isgn{E)sinZ ) ^' 

I 



where V-'s is the wavefunction for the scattering state, and bound states, respectively, p^"^ and 6p come from the 
Es and Eb are the energy of the scattering state and the scattering states, [s^] 
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The local density of states stemming from the bound 
states pb is given by 



Hamiltonian is 



Pb{E,x< {)) = - 



E ttmZ 
TTVp I cos Z I 



-2(_Btan Z)x/vf 



(9) 



for Eta.iiZ < and = for EtanZ > 0. The 
and 6p{E, x) are obtained in Ref. l32l while pb{E,x) has 
been missed in the literature. 
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FIG. 3: (Color online) (a) p{E,0)Ef, (b) p{E,-l/kF)EF, 
(c) pb{E,0)Ef and (d) pb{E, -l/kF)EF for several Z. Note 
Pb{E, x) — at Z = Q and also that the scales of (a), (c) and 
(b), (d) are different, and the contribution from ps is larger 
in (a) (at x — 0) compared with (b) (at x — — 1/A;_f). 



In Fig. H we show (a) p{E, 0)Ef, (b) p{E, -l/kF)EF, 
(c) pb{E,0)Ef and (d) pb{E, ~l/kF)EF for several Z. 
As shown in Fig. |3] (a) and Fig. [3] (b), the contribution 
from the bound states gives an asymmetry in the local 
density of states: pa{E,x) = p{E,x) — p{—E,x) ^ 
(note pb{E,x) = at Z = and also that p^^-' and 6p 
are even function of E). Since the bound states decay 
exponentially in space as seen in Eq. ^ , the contribution 
from the bound states becomes strongly suppressed and 
hence pa{E,x) goes to zero away from the gated region 
(compare Fig. [3] (c) and Fig. [3] (d)). The decay length of 
the bound states increases with the decrease of and 
therefore pB has a peak at low energy as shown in Fig. 
13] (d), which is reflected in the total density of states as 
shown in Fig. [3](b). Therefore, by comparing pa{E,x) 
near and away from the gated region using scanning tun- 
neling microscopy, one can identify the contribution from 
the bound states by Apa{E) = pa{E,0) — pa{E, — oo). 

Next, let us consider the tunneling current along the 
ID channel, taking into account the electron-electron in- 
teraction. Away from the half-filling, the possible scat- 
tering processes are dispersive (gd) and forward (g/) 
scatterings. 0, 01 Near the Fermi points, the correspond- 
ing interactions can be easily taken into account by 
the standard bosonization method. 3^37| The bosonized 



H^ - I dx 
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VF 



K = 



I 2ttvf |cos Z\+ gf — 2gd 
2ttvf |cos Z\ + gf + 2gd 



(10) 

(11) 
(12) 



with 
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^Li "^'"^ - ^flt ^ 'f'^i where (j)E.^ and 
(j)Li define chiral boson fields of spin up right mover and 
spin down left mover, respectively. As seen, the Luttinger 
parameter K is tunable by changing Z, namely the gate 
voltage. This opens up a possilbility of controlling the 
property of the ID interacting fermions by gate voltage. 

We consider a single magnetic impurity placed on the 
gated region, which produces the backward scattering. 
Then, the tunneling current Ic can be calculated within 
the linear response regime and is given 



Ic {eV)- 



(13) 



with an applied voltage V. Here, we consider a weak im- 
purity scattering. For a strong impurity scattering, the 
result will be modified due to formation of the localized 
impurity resonance. 28, Slj Thus, the dependence of K 
on Z can be experimentally studied through the tunnel- 
ing current. Since the backward scattering requires spin 
flip process, when the magnetic moment of the magnetic 
impurity is parallel to the spin of the bound states, the 
backward scattering is prohibited. This means that by 
changing the direction of the magnetic moment by, e.g. 
an applied magnetic field, the tunneling current can be 
tuned. 

In Fig. m we show (a) K and (b) 2/K — 1 as a function 
of Z/tt at gf/{2nvF) = gd/{2TTVF) = 0.1. As shown in 
Fig. |4](a), the K value oscillates with Z as expected from 
Eg. p^ . Near Z = tt/2 mod tt, K is strongly suppressed. 
Correspondingly, the exponent of the current, 2/K — 1, 
has a diverging behavior near such points as shown in 
Fig. m (b) . This dependence can be observed in tunneling 
experiments. 

It is noted that the conductance between the source 
and drain attached to the gated region is still dominated 
by the 2D surface metal on topological insulator, i.e., 
C2D = [e^ /h)kFl » cTiD ^ e'^/h with the mean free 
path £ and the conductance from the ID channel criu. 
One possible way to separate these two contributions is 
the modulation spectroscopy using the precession of the 
spin of the magnetic impurity. Namely, the backward 
scattering potential is dependent on the direction of the 
spin, which oscillates with the precession, and hence the 
ac component of the ID conductance occurs. One catch is 
that the universal conductance fluctuation (UCF) of the 
order oi e^/h also contributes to the ac component of the 
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FIG. 4: (a) K and (b) 2/;^ - 1 as a function of Zj-n at 
(7//(27ri;F) =pd/(27ruF) = 0.1. 



conductance [40|, but the dependence on the spin direc- 
tion is random, while the ID channel conductance shows 
the systematic dependence, i.e., it is most suppressed 
when the backward scattering is maximized. Therefore, 
by taking the average over several samples, one can sep- 
arate the UCF and the ID conductance of our interest. 

In summary, we studied the formation of the one- 
dimensional channels on the topological surface under 
the gate electrode. The energy dispersion of these chan- 
nels is almost linear in the momentum with the velocity 
sensitively depending on the strength of the gate voltage. 
The energy is also restricted to be positive or negative de- 
pending on the strength of the gate voltage. As a result, 
the local density of states near the gated region has an 
asymmetric structure with respect to zero energy. With 
the electron-electron interaction, the correlation effect in 
the bound states can be tuned by the gate voltage, which 
can be verified by tunneling experiments. 
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